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Abstract 

We investigate an exactly solvable model for directional transport in ID. The structured sys- 
tem, which has strong elastic interactions in its parts, explicitly demonstrates the role of volume 
exclusion in producing directional transport. We capture the complementary role of the elasticity 
and volume exclusion as the basic ingredient for showing up of broken microscopic symmetry at 
the scale of macroscopic motions. We compare the analytic results with the numerical simulation. 
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Brownian Ratchet (BR) paradigm finds wide applications when it comes to understanding 
the driven directional transport of microscopic objects There are a host of microscopic 

biological systems like kinesin and myosin motor proteins etc which move on polar tracks 
being driven by intermittent energy supply through ATP hydrolysis ^ . Attempts have been 
made to understand such motions in terms of BR paradigm invoking numerous variants of 
the BR model . Not only microscopic biological transport, there are many other areas 
like flow through micro channels, nanopores and even in cold atoms and optical lattices 
where the BR concept is applied to understand the transport phenomenon {4, ^. Brownian 

I I 

Motors - the term coined by Hanggi [4] for a class of such thermally assisted moving systems 
- require the system be driven out of equilibrium. This actually is a prerequisite of any 
directed transport. It also requires the symmetry be broken in some form either internally 
or externally. In many of these BR models, an internal degree of freedom of the microscopic 
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ect is coupled to an external field which is responsible for global symmetry breaking 



15| to result in macroscopic transport. In most of these models, the transported object 



is treated as a particle considering its internal structure not of importance. 

Recently, there has been an interesting report of a variant of BR model, so-called intrin- 
sic BR model, where the internal degree of freedom is not coupled to an external field for 
the sake of global symmetry breaking p^]. Rather, the internally broken symmetry at the 
microscopic scales shows up macroscopically in the transport of the object. The present 
work proposes an exactly solvable model that explicitly manifests reflection of microscopi- 
cally broken symmetry in the macroscopic transport of the system. We will be considering 
in the following a structured object as the motor which has symmetry broken inside and 
would not take into account any repeating ratchet potential (external field) that helps gen- 
erate global asymmetry. We will identify two important physical entities namely the elastic 
interaction between the parts of structured motor and volume exclusion as the essential 
and complementary requirements for the directed transport. While the elasticity allows the 
cross talk between the parts of the structured object, the volume exclusion is maintaining 
the locally broken symmetry intact and as a result of that the system shows macroscopic 
directional motion. We will also, substantiate our exact analytic results by comparing those 
with simulations. 

The present model takes into account, two particles connected by a spring and are im- 
mersed in a thermal bath characterized by a fixed temperature T. In such a situation the 
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model looks like 



xi = -a{xi - X2) + Tji 
X2 = —a{x2 — xi) + 772. 



(1) 



where, a is the force constant of the spring connecting the particles sitting at Xi and X2 
{xi > X2 say). The ?7's are the thermal noises experienced by the system where < r]i >= 
and < r}i{ti)r)j{t2) >— 2TSijd{ti — 12) which is the fluctuation-dissipation relation. Note 
that, the a and the Boltzmann constant kb have been absorbed by adjusting the scale of 
the temperature in the fluctuation-dissipation relation which is a very standard practice. 

Such a system would soon attain equilibrium and one naturally cannot expect having di- 
rectional transport. So, wc put a drive to one of the particles which will be driving the system 
out of equilibrium by breaking the fluctuation-dissipation relation. Throughout this paper, 
we will consider an over damped dynamics of the system since this over damped regime 
is physically relevant in most of the biological and other cases. Now, in the overdamped 
regime, where there is no inertial term present, the velocity of a particle is proportional to 
the force acting on it. Taking this into consideration we will use a drive which is equivalent 
to feeding back the instantaneous velocity to one of the particles and that would invalidate 
the otherwise existing fluctuation dissipation relation. This not only is a very general way of 
driving the system out of equilibrium, but also it would enable us to solve the system exactly 
and calculate the average velocity and efficiency of the system in exact analytic forms. So, 
with this driving, our model looks like 



where /3 is the velocity feedback constant or the strength of the external drive. In what 
follows we will keep — 1 < /3 < 1 in order to remain safely within the over damped regime. 
A large forcing would eventually require the inertial term be taken into consideration which 
we would not like to move into at the moment. 

In order to understand where this modified situation with the velocity feedback actually 
differs from the model in Eq.l let us write it down in an analogous from as 



xi = -a{xi - X2) + /3xi + rji 
X2 = —a{x2 — xi) + 772. 



(2) 



X2 = -Ol(x2 - Xi) + 772. 
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(3) 
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the velocity feedback term which we are considering here as a drive (which one might as well 
call the extra damping on the particle at xi than that in the position X2) is in effect changing 
the local temperature of the particle at Xi than that at X2 resulting in a nonequilibrium 
situation. This nonequilibrium situation would sustain all over the time because within the 
scope of the model the Gaussian noise strength and the so called damping on the particles are 
held constant. Thus, effectively in this model the local environment of the two interacting 
particles are kept fixed at two different temperatures and that keeps the nonequilibrium 
situation prevailing for all the time. The above reasoning makes it clear why, in the present 
context, we prefer calling the feedback term as a drive than considering it merely a part of 
the damping on the particle at Xi. 

To analyze this model further, let us move to the center of mass (CM) X = {xi + X2)l2 
and the internal gap Z = {x\— X2) coordinates. The model now reduces to the form 

. a(2-/3) 



1-/3 



z 



where, < >= 0, < Czmz{t') >= 2TD^5{t - t'), < Cxmx{t') >= ^S{t - t') and 

Eq.2 clearly shows that the internal coordinate Z will have an equilibrium distribution 

P^AZ) = expi ^^^^=^^^2). (5) 

Since, the average velocity of the CM < X > is proportional to the < Z >, it can be 
easily calculated. Note that, although Z has a probability distribution with zero mean the 
volume exclusion will prevent Z being less than zero and thus we will actually have a nonzero 
< Z > by simply integrating one one half of the distribution it has. This is exactly where 
we explicitly see the role of volume exclusion in producing directional transport. Thus, we 
arrive at 



<Z> 



l2TD,il-(3) 
7ra(2 - /3) 



^'''^-"^(W)^^^-"^(W)"V -(2-/5) • 
Note that, at ^ = the average velocity of the system vanishes identically. This is in keeping 
with the fact that one cannot expect directed transport in equilibrium. Interestingly, in the 



above expression, we see that with a positive ^ the system will move towards left i.e. in 
the direction from xi to X2 and it changes sign with the change in sign of ^. This selection 
of the direction of motion can easily be understood. The effective spring constant for the 
particle at x-y being a/(l — /3) is greater than a for a positive /3. As a result the particle at 
Xi will impinge on the particle at X2 with more momentum and that will make the system 
move towards left. The situation will change with the change in the sign of j3. 

A calculation of the efficiency of the system is quite straight forward. Considering the 
average input energy to the system is proportional to < X"^ > (i.e. average instantaneous 
kinetic energy of the CM), we can take the effective energy output as < AT >^. Now, the 
efficiency can be given as rjeff =< AT >^ / < >. The explicit calculation produces 

Some interesting features can immediately be found out from the above relation for the 
efficiency. The efficiency being not dependent on the temperature it explicitly shows that 
the thermal energy is not being drained out. But, it is helping the system achieve directional 
transport because, the average velocity of the CM is a function of temperature. This explicit 
non-dependence of the efficiency on temperature is an important feature. The maximum 
efficiency (for ^ = 1) is 100 x 2/7r which is roughly 64%. This is a big percentage efficiency 
indeed. The efficiency goes linearly to zero with a and with (3. 

To compare the analytic results we simulate the system. While simulating the system, to 
apply the volume exclusion i.e. correcting for particle positions when they cross each other 
during sequential updating we work at two extreme limits namely particles suffering totally 
elastic and inelastic collisions. In the totally elastic collision scheme, we swap the particle 
positions following any violation of the volume exclusion and thus the CM remains unaltered 
preserving the momentum. In the totally inelastic scheme, we place the particles at the CM 
following any violation of the volume exclusion constraint. In Fig. la and b we have plotted 
the average velocity of the system against a and /3 respectively at some constant value of 
the other parameters. During these simulations we have kept the strength of the thermal 
noise < 77^ >= 2T — 0.0001. In all the simulations we have taken the maximum positive 
value for /3 to be 0.99 because the average velocity of the system diverges at ^ = 1. In 
Fig. 2a and b we have plotted the efficiency of the system against the parameters a and P 
as well under similar situations. In these figures a comparison of simulation results with 
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theoretically calculated numbers is excellent and directly justifies our exact calculations. 

In most of the BR models actual structural realities like volume exclusion etc are masked 
due to considering the system as a structureless particle. Moreover, considering a repeated 
setup of asymmetric microscopic potentials in ID is also in effect a macroscopic symmetry 
breaking that results in macroscopic directional transport. With the help of the present 
model, we explicitly show that in ID there is no need for such macroscopically broken 
symmetry present for the transport to happen. Considering the internal structure of the 
microscopic moving system, we show that the volume exclusion would actually make the 
microscopic broken symmetry remain intact over all the time in ID where there is no ro- 
tational degrees of freedom. Thus, the direction of motion picked up by the system will 
be there for ever if the nature of the instantaneous velocity feedback remains the same. In 
this connection a discussion of the role of the feedback term in our model is worth doing. 
In effect, the feedback term is serving two purposes - it breaks the symmetry and drives 
the system out of equilibrium. Actually, for /3 > 0, 1. it reduces the damping experienced 
by the particle at xi and 2. it makes the temperature of the particle at xi different from 
that at X2 (which is a non-equilibrium situation ensuring possibility of heat flow between 
the parts of the system). A change in the sign of /3 would do the same but in a way so that 
the particle at X2 would now be in the driver's seat by having smaller damping and larger 
temperature. The different damping experienced by two different particles could be because 
of two differently shaped bodies sitting at xi and X2 and the difference in temperature could 
be due to the fact that one of them is being driven by some unbiased force. Within the 
present feedback picture, we do not need to take the shape of the objects at positions xi 
and X2 into consideration, rather take them as particles. However, taking these two essen- 
tial requirements like symmetry breaking and the driving in a correlated way through the 
velocity feedback is obviously a restriction. But, it nicely allows us to do a simple exact 
analytic calculation for the system to explicitly understand the interplay of the elasticity 
and volume exclusion to generate directional transport. An extension of the system in 2D, 
where rotational degrees of freedom are present, is straight forward 171] and it can generally 
capture the principle of transport of some bacterial systems. 
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FIG. 1: a. Average velocity of the CM vs a. 
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FIG. 1: b. Average velocity of the CM vs j5. 
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FIG. 2: a. Efficiency [r]) vs a. 
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